Dipolar Bose-Einstein Condensates with Weak Disorder 
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A homogeneous polarized dipolar Bose-Einstein condensate is considered in the presence of weak 
quenched disorder within mean-field theory at zero temperature. By first solving perturbatively the 
underlying Gross-Pitaevskii equation and then performing disorder ensemble averages for physical 
observables, it is shown that the anisotropy of the two-particle interaction is passed on to both the 
superfluid density and the sound velocity. 
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Bose-Einstein condensates (BECs) in a disordered en- 
vironment have been the subject of various experimental 
investigations in recent years. Historically, this "dirty 
boson" problem arose in the context of superfluid helium 
in Vycor glass Later, disorder either appeared nat- 
urally, for instance in magnetic wire traps [2|, [3] , or was 
created artificially and controllably by using laser speck- 
les 0, [B| . Most theoretical studies focused on analyzing 
the Bogoliubov theory of dirty bosons for weak disorder 
within either the second quantization or the replica 
method , [Til ] . In this way, it turned out that superflu- 
idity persists despite quenched randomness, but a deple- 
tion occurs due to the localization of tiny condensates in 
the respective minima of the disorder potential. So far, 
this localization scenario has only been analyzed for the 
contact interaction, as it is usually dominant for ultra- 
cold dilute Bo se g ases. Since the realization of atomic 
dipolar BECs [12j and the generation of heteronuclear 
molecules in the rovibrational ground state near quan- 
tum degeneracy long-range and anisotropic dipole- 
dipole interactions have also attracted attention [l4-ll6|. 
Therefore, we consider in this Rapid Communication the 
impact of weak disorder upon a polarized dipolar BEC at 
zero temperature. In particular, we will show that both 
the superfluid density and the sound velocity yield char- 
acteristic interaction-induced anisotropies, which are not 
present at zero temperature in the absence of disorder. 

The weakly interacting theory of dirty bosons 
11] states that Bogoliubov quasiparticles and disorder- 
induced fluctuations decouple in the lowest order. This 
suggests a simplified approach in which the leading cor- 
rection due to the presence of a random potential could 
be derivable from a mean-field theory. Therefore, we as- 
sume at T = that all bosons occupy the same quan- 
tum state, for which the macroscopic wave function ^(x) 
obeys the time-independent Gross-Pitaevskii (GP) equa- 
tion 
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U( X )+ J dV|*(x')|V int (x-xO 
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and neglect with this from now on any impact of quantum 
fluctuations. Here, m denotes the particle mass, /i stands 
for the chemical potential, V; n t(x— x') represents an arbi- 
trary two-particle potential with inversion symmetry, and 
U(x) describes the disorder potential, which is defined by 
its statistical properties. Denoting the disorder ensem- 
ble average according to ( • ) , a homogeneous disordered 
system has a vanishing first moment (C/(x)) = and the 
second moment is of the form ([/(x)[/(x')) = R(x — x'). 
The GP equation (fTJ) represents a stochastic nonlinear 
partial differential equation, where the given statistics 
of the disorder potential t/(x) are mapped to the wave 
function \P (x) [17| . In the case that the random potential 
U(x) is small in comparison with all other energy scales, 
its perturbative treatment is justified. To this end, we 
decompose the wave function of the system according to 



*(x) = Vo(x) + Vi(x) + $j(x) 



(2) 



and solve the GP equation in the zeroth, first, and second 
order of U(x), respectively. As the expansion (|2J) deter- 
mines the ground state, the wave function Vl'(x) turns out 
to be real. Afterwards, we determine the disorder ensem- 
ble average for both the particle density n = {^(x) 2 ) and 

the condensate density no = (^(x)) 2 , thus, the conden- 
sate depletion results in the lowest order in 



n 



d 3 k 



R(k) 



(2tt)3 [/j2 k 2/2m + 2nT/ int (k)] 
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Note that the existence of the k integral implies that the 
Fourier transform of the interaction potential V{ n t (k) has 
to be strictly positive. Physically, this condensate deple- 
tion is due to the formation of fragmented condensates in 
the respective minima of the random potential. In order 
to further quantify this notion, a separate order param- 
eter was recently proposed in Ref. [18|. It is motivated 
by defining the condensate density no as usual as the su- 
perfluid order parameter from the off-dia gon al long-range 
order of the one-particle density matrix [19( 
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The density of fragmented condensates q is then identi- 
fied as a separate Bose-glass order parameter, similar to 
the Edwards- Anderson order parameter of a spin glass 
20]. It can be defined by considering the off-diagonal 
long-range order of the two-particle density matrix 



lim (*(x) 2 *(x') 2 ) = {n Q +q) 2 , 

|x— x'|— >oo 
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as the latter has to coincide with the square of the total 
density n. Applying the concept of Eqs. (HJ) and §5§ to the 
perturbative solution ((2|) of the GP equation ([1]) indeed 
yields, together with Eq. ([3]), the result that the density 
of the fragmented condensates q defined in Eq. ([5]) coin- 
cides in the lowest order with the condensate depletion 
ri — n found in Eq. Thus, we conclude that the 

localization phenomenon for weak quenched disorder fol- 
lows already from a mean-field description of the dirty 
boson problem. Therefore, our mean-field approach rep- 
resents a simplified derivation for the disorder-induced 
condensate depletion Q in comparison with the Bogoli- 
ubov theory of Refs. [q-lllj. Note that, recently, disorder 
effects for Bogoliubov quasiparticles have been analyzed 
in Ref. [p. 

For a polarized dipolar BEC, the particles interact via 
a contact interaction with strength g — 4irh 2 a/m, where 
a denotes the s-wave scattering length, and they pos- 
sess dipole moments, which are aligned along the z axis. 
Thus, the Fourier transform of the interaction potential 
is given by Mnt(k) = g[l + edd(3a: 2 — 1)], where we define 
x = k • e z = cosr? and edd = Cdd/3g with Cdd = Mo m2 , 
in the case of magnetic dipoles with the dipole moment 
m, and Cdd = 47rd 2 for electric dipoles, where the dipole 
moment d is measured in Debyes. In the case of a spa- 
tially decaying disorder correlation i?(x), the results of 
the theory do not depend significantly on its shape 22 1. 
Therefore, in what follows, we restrict ourselves to the 
case of a Gaussian correlation with the Fourier trans- 
form i?(k) = Re~ a k / 2 , where R and a characterize 
the strength and the correlation length of the disorder, 
respectively. With this, the condensate depletion (|3|) spe- 
cializes to ti — hq — tihm / (edd ) + • • • j where the limit 
of pure contact interaction and delta-correlated disorder 
yields n RM = [m 2 R/(8TT 3 ^h 4 )} ^fnfa 0. Whereas the 
relative interaction strength edd increases the condensate 
depletion, the ratio of the correlation length a and the 
coherence length £ = 1/ V8vma decreases it according to 
the function 
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[I + 2C(x)]e c ^VfcVC(20 2cr 

\fiE dd X 2 + 1 - £ d d V^C 
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with the abbreviation ((x) = <r 2 (3edd^ 2 + 1 — £dd)/£ 2 
(see Fig. [1]). In particular, for small £dd, we have 
/ (edd, <r/0 = A{a/0 + B(a/0 e 2 d + . . ., with positive 
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FIG. 1: (Color online) Dipolar condensate depletion func- 
tion ([6]) with the special cases of vanishing correlation length 
/(edd,0) and pure contact interaction /(0, <x/£) H indicated 
by solid lines. 



coefficients A(cr/£) and B(a/£) which decrease with in- 
creasing disorder correlation length. In contrast, the con- 
densate depletion function f (edd, c/f) diverges in the 
limit edd t 1) i n accordance with the conclusion that 
the existence of the k integral in Eq. ((3]) is lost provided 
Vint (k) vanishes for any k. This indicates that our pertur- 
bative treatment breaks down provided that the dipolar 
interaction is strong enough relative to the contact inter- 
action. 

In order to describe a superfluid within mean-field the- 
ory, we apply a Galilean boost with wave vector kN to 
the time-dependent GP equation and introduce a moving 
condensate via the ansatz Vl'(x) = ip(x.)e lks ' x (23j |. With 
this, we obtain in the stationary case 



h 2 H 2 
— A - i — K • V + U(x) - fi eS 
2m m 



+ / dVVW(x-x')|V(x')| 2 



V(x) = 0, (7) 



where we have introduced the abbreviations K = ks — kfj 
and fi c ff = /i — S- 2 k|/2m + h 2 h^ ■ ks/rn. Within a linear 
response, we can assume kN and ks to be small, so that 
the superfluid part of the system will move with a wave 
vector proportional to ks and the nonsuperfluid part is 
coupled to k]^. Therefore, the expansion of the total 
momentum P = — J rf 3 i$*(x)iSV$(x) with respect to 
ks and kw yields in the disorder-averaged case (P) = 
V (nsfiks + riN^-kN) + . . . , where V denotes the volume 
of the system and n§ and «n represent the superfluid 
and the nonsuperfluid density, respectively. Expanding 
the solution ip(x) of Eq. (JTJ), which may now be complex, 
with respect to the disorder yields in general the result 
that the superfluid density turns out to be a tensor with 
the components 



d 3 k 



4nR(k)kikj 



(2tt) 3 k 2 [h 2 k 2 /2m + 2nV int (k)] 



.(8) 



and ns + tin = n. In the case of isotropy, i.e. i?(k) = 
i?(|k|) and T4nt(k) = V[nt(|k|), this reduces to a diagonal 
superfluid density ns — n — 4(n — no)/ '3 + . . „ which 
generalizes the case of pure contact interaction [6l-lll [ . 



3 




FIG. 2: (Color online) Superfluid depletion function a) parallel and b) perpendicular to the polarization axis of the dipoles. 
The solid lines indicate the special cases of vanishing correlation length and pure contact interaction, respectively. The dashed 
line defines the critical relative interaction strength via / (ead.crit, °7£) — /|| ( e dd,crit, cr/£). 



For a dipolar BEC and Gaussian disorder correlation 
function, Eq. (J5|) yields a superfluid density that depends 
on the direction of the superfluid motion with respect to 
the orientation of the dipoles. Parallel to e z , the super- 
fluid depletion reads n — n s n = tihm/|| (edd)C/£) + . . . 
with the function 



f\ 



£dd, j ) =4 



dx 



x 2 [l + 2C{x)]e^eric^/0x) 8c 



\JieadX 2 + 1 - £dd 



3y^ 
u 

(9) 

whereas, perpendicular to the dipoles, we have n— ns,± — 
«hm/i (fdd, <t/Q + ■■■ with the function f± (e dd , a/0 = 
2/(edd,o"/0 - /|| (fdd,cr/0 /2. For small e d d, we get 
/||,_L (edd,cr/0 = % x (ct/0 + £ dd + ■ • ■ , again 

with coefficients A\\_^(a/0 and B\\_^{a/Oi whose abso- 
lute values decrease with increasing disorder correlation 
length. This time, however, only f± (e dd ,<r/£) diverges 
for e dd t 1) whereas /j| (e dd ,er/£) remains finite in this 
limit. 

Comparing Figs. [2] a) and[2]b) by taking into account 
the different vertical scales, we conclude that the super- 
fluid density parallel to e z is always less depleted than 
the superfluid density perpendicular to e 2 . This result 
can be qualitatively explained by considering Vj nt (k), in- 
cluding the dipolar interaction, as an effective contact 
interaction strength for a flow with wave vector k. This 
quantity is smaller than g perpendicular to e z and larger 
than g parallel to e z . As the Huang-Meng depletion uhm 
for a pure contact interaction scales with the power 
— 1/2 with respect to the contact interaction strength, 
we obtain in this picture a smaller depletion parallel 
to e z than perpendicular to e z . In addition, we note 
that /|| (0, a/0 > f (0, a/0 |, whereas (e dd , a/0 de- 
creases and /(e dd ,cr/£) increases, respectively, with e dd 
for fixed a/0 Therefore, a system with a sufficiently large 
relative interaction strength e dd has the astonishing prop- 
erty that the depletion of the parallel superfluid density 
is smaller than the condensate depletion. Figure [2] a) re- 
veals in the dashed line how the critical value e dd-cr i t de- 
creases with increasing disorder correlation length a. The 
existence of e ddjCr it represents a counterintuitive result, as 
particles of the fragmented BECs, which are supposed to 



be localized in the respective minima of the random po- 
tential, seem to contribute to the parallel superfluid mo- 
tion of the system. We conclude that, due to the presence 
of the dipolar interaction, the locally condensed particles 
are only localized for a certain time scale. For longer 
time periods, our finding suggests that an exchange of 
the localized particles occurs with the nonlocalized par- 
ticles, thus allowing for a superfluid density that is larger 
than the condensate density. This supports the finding 
of Ref. 18(, in which such a finite localization time for 
tiny BECs was calculated within a Hartree-Fock theory 
of dirty bosons for the special case of a pure contact in- 
teraction. 

Finally, we determine the speed of sound within a hydro- 
dynamic approach [24|. To this end we assume that, in 
the long-wavelength limit, the respective transport quan- 
tities are effectively disorder averaged. Thus, we consider 
both the Euler equation mdvs/dt + V Ut + mvg/2) = 
and the continuity equation dn/dt + V • j = with the 
current density j = ns v s + "nvn, where vs and vn de- 
note the superfluid and boost velocity, respectively. As 
the normal component is pinned, we must also assume 
that it remains stationary. This leads to the condition 
vn = 0, which is reminiscent of the physically closely 
related problem of fourth sound in 4 He [24j . Considering 
small space- and time-dependent perturbations around 
the equilibrium values vg — <5vg(x, t), n — n cq + Sn(x, t), 
n s = ns,eq + fins (x, t), and /i = jj,[n eq + Sn(x, t)] yields 



d 2 6n{x,t) 
dt 2 



n s , cq dp, 
m dn 



V5n(x,t) 



= 0. 



(10) 



With this, we obtain in the sound wave regime for the 
speed of sound in the direction q the general result 
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d 3 k 



i?(k) 



^nV int (0)/m ' ' ~J (2^) 3 [h 2 k 2 /2m + 2nV in t(k)} 2 
H 2 k 2 V int (k)/2mV int (0 



h 2 k 2 /2m + 2nV int (k) 



- mr 



(ii) 



where the first term originates from the equation of state 
and the second term originates from the superfluid den- 
sity ©. Note that our hydrodynamic derivation of the 
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FIG. 3: (Color online) a) Speed of sound for vanishing disorder and b) correction term s(tdd,'&) for nonvanishing disorder in 
Eq. (E). 



speed of sound reduces for the special case of a contact 
interaction to an expression that has recently been ob- 
tained within an independent Green's function approach 
25j . In the case of a dipolar interaction and an uncorre- 
cted disorder i?(x) = RS(x), we have 



c(£dd,^) 
co(edd,tf) 



1 



2c2(e dd ,i?) 



(edd,0) 



(12) 



with co(edd,$) = y/ngfadd cos 2 d + 1 — eddj/m repre- 
senting the speed of sound for vanishing disorder and 
with a dipolar function s(edd,$), which generalizes the 
result for the pure contact interaction s(0, #) = 5/3 0- 
Figure [3] shows that due to weak disorder, the speed of 
sound has a strong dependence on both the sound di- 
rection, characterized by q • e z — cosi9, and the relative 
interaction strength tdd- 

Finally, we conclude that the delicate interplay of dipo- 
lar interaction and weak disorder yields characteristic 
anisotropies for physical observables at T = 0. The 
anisotropy of the speed of sound should be detectable 
with modern Bragg spectroscopy by measuring the un- 
derlying dynamic structure factor 26|, |27| . Even more 
important could be the anisotropy for the superfluidity 
density, as this should affect the collective excitations 
[28j for harmonically trapped dipolar condensates within 
a random environment. It is expected that these inter- 
esting effects are more pronounced for strong dipolar in- 
teractions that arise for highly magnetic atoms, such as 
dysprosium 



2911 . or polar heteronuclear molecules, such 
13J. To this end, it is indispensable to in- 



as 40 K 87 Rb 

vestigate dirty dipolar BECs in confined geometries, for 
instance in a harmonic trap 30] or an optical lattice 31 1 . 

We thank J. Dietel, R. Graham, M. von Hase, 
P. Navez, and F. Nogueira for useful discussions and ac- 
knowledge financial support from the German National 
Scholarship and from the German Research Founda- 
tion (DFG) within the Collaborative Research Center 
SFB/TR12. 



[2 
[3 

[4 

[5. 
[6 

[T 

[8 

[9 

[io; 
in 

[12 

[13 
[14 
[15 
[16 

[17; 

[18 
[19 

[20 

[21 

[22 

[23 

[24 
[25 
[26 



[1] B. C. Crooker, B. Hebral, E. N. Smith, Y. Takano, and 
J. D. Reppy, Phys. Rev. Lett. 51, 666 (1983). 



P. Kriiger et al, Phys. Rev. A 76, 063621 (2007). 

J. Fortagh and C. Zimmermann, Rev. Mod. Phys. 79, 

235 (2007). 

J. C. Dainty (editor), Laser Speckle and Related Phenom- 
ena (Springer, Berlin, 1975). 

J. Billy et al, Nature (London) 453, 891 (2008). 

K. Huang and H.-F. Meng, Phys. Rev. Lett. 69, 644 

(1992). 

S. Giorgini, L. Pitaevskii, and S. Stringari, Phys. Rev. B 
49, 12938 (1994). 

M. Kobayashi and M. Tsubota, Phys. Rev. B 66, 174516 
(2002). 

G. M. Falco, A. Pelster, and R. Graham, Phys. Rev. A 

75, 063619 (2007). 

A. V. Lopatin and V. M. Vinokur, Phys. Rev. Lett. 88, 
235503 (2002). 

R. Graham and A. Pelster, in Path Integrals - New 
Trends and Perspectives: Proceedings of the 9th Inter- 
national Conference, edited by W. Janke and A. Pelster 
(World Scientific, Singapore, 2008), p. 376. 
A. Griesmaier, J. Werner, S. Hensler, J. Stuhler, and T. 
Pfau, Phys. Rev. Lett. 94, 160401 (2005). 
K. K. Ni et al, Science 322, 231 (2008). 
M. A. Baranov, Phys. Rep. 464, 1 (2008). 
L. D. Carr and J. Ye, New J. Phys. 11, 055009 (2009). 
T. Lahaye, C. Menotti, L. Santos, M. Lewenstein, and T. 
Pfau, Rep. Prog. Phys. 72, 126401 (2009). 
P. Navez, A. Pelster, and R. Graham, Appl. Phys. B 86, 
395 (2007). 

R. Graham and A. Pelster, Int. J. Bif. Chaos 19, 2745 
(2009). 

A. J. Leggett, Quantum Liquids: Bose Condensation and 
Cooper Pairing in Condensed Matter, (Oxford University 
Press, New York, 2006). 

S. F. Edwards and P. W. Anderson, J. Phys. F 5, 965 

(1975). 

C. Gaul and C. A. Miiller, e-print larXiv:1009. 54481 

M. Timmer, A. Pelster, and R. Graham, Europhys. Lett. 

76, 760 (2006). 

M. Ueda, Fundamentals and New Frontiers of Bose- 
Einstein Condensation (World Scientific, Singapore, 
2010). 

I. M. Khalatnikov, An Introduction to the Theory of Su- 
perfluidity (Benjamin, New York, 1965). 
C. Gaul and C. A. Miiller, Phys. Rev. A 83, 063629 
(2011) 

J. Stenger, S. Inouye, A. P. Chikkatur, D. M. Stamper- 
Kurn, D. E. Pritchard, and W. Ketterle, Phys. Rev. Lett. 
82, 4569 (1999). 



■5 



[27] R. Ozeri, N. Katz, J. Steinhauer, and N. Davidson, Rev. 

Mod. Phys. 77, 187 (2005). 
[28] G. M. Falco, A. Pelster, and R. Graham, Phys. Rev. A 

76, 013624 (2007). 
[29] M. Lu, S. H. Youn, and B. L. Lev, Phys. Rev. Lett. 104, 



063001 (2010). 
[30] G. M. Falco, J. Phys. B 42, 215303 (2009). 
[31] P. Muruganandam, R. K. Kumar, and S. K. Adhikari, J. 

Phys. B 43, 205305 (2010). 



